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ABSTRACT 
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O \ We present photometry of HD 189733 during eight transits of its close-in giant planet, 
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and out-of-transit photometry spanning two years. Using the transit photometry, we 
determine the stellar and planetary radii and the photometric ephemeris. Outside of 
{Sj ■ transits, there are quasiperiodic flux variations with a 13.4 day period that we attribute 

to stellar rotation. In combination with previous results, we derive upper limits on the 
j> ■ orbital eccentricity, and on the true angle between the stellar rotation axis and planetary 

orbit (as opposed to the angle between the projections of those axes on the sky). 
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1. Introduction 

o ■ 

For the same reason that eclipsing binary stars are important in stellar astrophysics, transiting 
planets play an outsized role in exoplanetary science. This can be appreciated by comparing the 
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transiting exoplanet HD 189733b (Bouchy et al. 2005), the subject of this paper, to the arbitrarily- 
chosen non-transiting planet HD 187123b (Butler et al. 1998). Both planets are "hot Jupiters" 
detected by the Doppler method. All that is known about HD 187123b is its orbital period (P = 
3.097 days) and minimum mass (M p sini = 0.52 Mj up ), despite 6 years having elapsed since its 
discovery. In contrast, for HD 189733b, discovered only 1.5 yr ago, transit photometry has revealed 
the planet's radius (1.15 i?j up ; Bouchy et al. 2005, Bakos et al. 2006a) and removed the sini 
ambiguity in the planet's mass (M p = 1.13 Mj up ). Infrared photometry during a secondary eclipse 
has led to a determination of the planet's 16 /im brightness temperature (T p = 1300 K; Deming 
et al. 2006). Most recently, spectroscopic observations during a transit have shown that the angle 
on the sky between the orbit normal and the stellar rotation axis is within a few degrees of zero 
(Winn et al. 2006). 

These measurements are essential for a complete understanding of the atmospheres and inte- 
riors of hot Jupiters, as well as their formation and migration mechanisms. An even richer set of 
measurements can be expected in the future, from investigators pursuing transmission spectroscopy, 
reflected- light observations, and other transit-related investigations (as reviewed recently by Char- 
bonneau et al. 2006a). One goal of the Transit Light Curve (TLC) Project is to support these 
efforts by refining the estimates of the planetary, stellar, and orbital parameters, through high- 
accuracy, high-cadence photometry of exoplanetary transits. We also seek to measure or bound 
any variations in the transit times and light-curve shapes that would be caused by the influence of 
additional bodies in the system (Miralda-Escude 2002; Agol et al. 2005; Holman & Murray 2005). 
Along the way, we are exploring different techniques for photometry and parameter determina- 
tion. Previous papers in this series have reported results for the exoplanets XO-lb (Holman et 
al. 2006a), OGLE-TR-lllb (Winn et al. 2007), TrES-1 (Winn, Holman, & Roussanova 2006), and 
OGLE-TR-lOb (Holman et al. 2006b). 

This paper presents our results for HD 189733b, along with out-of-transit photometry spanning 
two years. The reason for gathering out-of-transit photometry was to attempt to measure the stellar 
rotation period. The parent star is relatively active, with a chromospheric activity index S = 0.525 
(Wright et al. 2004) and log R' HK = —4.4, raising the possibility of measuring the rotation period 
through starspot-induced quasiperiodic flux variations. As we will explain, the measurements of 
both the rotation period and the Rossiter-McLaughlin effect permit the determination of the true 
angle between the orbit normal and stellar rotation axis (as opposed to the angle between the 
projections of those vectors on the sky, which was measured by Winn et al. 2006). 

This paper is organized as follows. In § 2, we present photometry of 8 different transits, along 
with nightly photometry over two consecutive observing seasons. In particular, § 2.2 presents the 
out-of-transit photometry and the estimation of the stellar rotation period. In § 3, we describe the 
parameteric model that was fitted to the data, and in § 4 we present the results for the planetary, 
stellar, and orbital parameters, as well as the new transit ephemerides, a limit on the orbital 
eccentricity, and the three-dimensional spin-orbit alignment. The last section is a summary. 
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2. Observations and Data Reduction 

Our observations took place in 2005 and 2006 using telescopes at 4 different observatories. 
Each of the sections below presents the photometry from a given observatory. All together, we 
observed 8 different transits. Table 1 gives a summary of the characteristics of the data from each 
transit, and Table 2 gives the final photometry. A telescope at Fairborn Observatory (§ 2.2) was 
also used to monitor the out-of-transit flux of HD 189733 over the two observing seasons. 

2.1. Fred Lawrence Whipple Observatory 

We observed the transits of UT 2006 Jul 21 and Sep 10 with the 1.2m telescope at the Fred 
L. Whipple Observatory on Mt. Hopkins, Arizonafj] We used KeplerCam, which has one 4096 2 
Fairchild 486 back-illuminated CCD, with a 23'1 x 23'.1 field of view. For our observations we used 
2x2 binning, which gives a scale of 0768 per binned pixel, a readout and setup time of 11 s, and 
a typical readout noise of 7 e~ per binned pixel. We observed both transits through the SDSS 
z band filter, in order to minimize the effect of color-dependent atmospheric extinction on the 
relative photometry and to minimize the effect of limb-darkening on the transit light curve. We 
deliberately defocused the telescope such that the full-width at half-maximum (FWHM) of stellar 
images was about 7 binned pixels (4'/8) in order to permit a consistent exposure time of 5 s. We 
used automatic guiding to keep the image registration as constant as possible. We also obtained 
dome flat exposures and zero-second (bias) exposures at the beginning and the end of each night. 

On UT 2006 Jul 21, the sky conditions began partly cloudy but gradually improved as the 
night went on. We observed the target star as it rose from an airmass of 1.19 to 1.01 and then 
descended to an airmass of 1.23. There was a 15-minute interruption after second contact, due to 
clouds. On UT 2006 Sep 10, the sky conditions were mainly clear (but not all-sky photometric). 
We followed the target star from an airmass of 1.02 to 2.5, although we discarded the data taken 
at airmass > 1.9 because of their much poorer quality. There was a 15-minute interruption prior 
to third contact, due to a computer crash, which also caused a change in image registration and 
focus. 

We used standard IRAfU procedures for the overscan correction, trimming, bias subtraction, 
and flat-field division. We performed aperture photometry of HD 189733 and 14 nearby and 
necessarily fainter stars. The light curve of each comparison star was normalized to have unit 
median, and the mean of these normalized light curves was taken to be the comparison signal. 
The light curve of HD 189733 was divided by the comparison signal, and corrected for residual 



lr The data from the first of these two transits have already been presented by Winn et al. (2006). 

2 The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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systematic effects by dividing out a linear function of time. The zero point and slope of the linear 
function were determined as part of the model-fitting procedure, as explained in § 4. Figure 1 
shows the final light curves, along with a time-averaged composite light curve created from the two 
data sets. 

2.2. Fairborn Observatory 

We used the T10 0.8m automated photometric telescope (APT) at Fairborn Observatory, in 
Arizona, to observe four complete transits of HD 189733b and to monitor the out-of-transit stellar 
flux. The T10 APT is equipped with a two-channel precision photometer employing two EMI 
9124QB bi-alkali photomultiplier tubes to make simultaneous measurements in the Stromgren b 
and y passbands. The APT measures the difference in brightness between a program star and 
a nearby constant comparison star (or stars) with a typical precision of 0.0015 mag for bright 
stars (V < 8.0). For the HD 189733 transits, we used the comparison star HD 189410 (V = 5.68, 
B — V = 0.34, F0). The differential magnitudes were reduced with nightly extinction coefficients and 
transformed to the Stromgren system with yearly mean transformation coefficients. To improve our 
photometric precision, we combined the separate b and y differential magnitudes into a single (b + 
y)/2 passband. For additional information on the telescopes, photometers, observing procedures, 
and data reduction techniques, see Henry (1999) and Eaton, Henry, & Fekel (2003). 

We observed the transits of UT 2005 Nov 28, 2006 May 2, 2006 May 22, and 2006 Jun 11. 
On each of those nights, the differential magnitude of HD 189733 was recorded for 3-5 hours 
bracketing the expected time of mid-transit. The transit light curves are shown in the top 4 panels 
of Fig. 2. Like the FLWO data, these data have been corrected by a linear function of time that 
was determined as part of the fitting procedure (see § 4) . 

In addition, we also measured the out-of-transit flux with the APT on 93 different nights 
spanning two observing seasons between 2005 October and 2006 July. Three comparison stars were 
observed on each night. The out-of-transit flux measurements are shown in the top two panels 
of Fig. 3. The flux varied erratically during the 2005 observing season and at the beginning of 
the 2006 season. However, a quasiperiodic signal became evident near the end of our 2006 season 
observations with a peak-to-peak amplitude of 1.3% and a period of about 13 days, although only 
2.5 cycles were observed. 

This type of photometric behavior — erratic and occasionally quasiperiodic variation — is com- 
mon among chromospherically active stars, especially stars with low or intermediate levels of ac- 
tivity (see, e.g., Henry, Fekel, & Hall 1995). The photometric variability in these stars arises from 
photospheric starspots and plages that are carried into and out of view by the stellar rotation. In 
the case of HD 189733, the spots cover only ~1% of the stellar surface at any moment, and their 
distribution on the star changes significantly on the rotational timescale. A periodogram of the 
quasiperiodic portion of the light curve (indicated by the solid dots in the middle panel of Fig. 3) 
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Fig. 1. — Photometry of HD 189733 in the z band, using the FLWO 1.2m telescope and Keplercam. 
These data were used to estimate the planetary, stellar, and orbital parameters (see § 3). The 
bottom panel is a composite light curve created from both data sets, after time-shifting the earlier 
transit and averaging into 1 minute bins. The residuals (observed— calculated) are plotted beneath 
the data. 
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shows a strong peak at 13.4 days, which we take to be the stellar rotation period. We estimated 
the uncertainty by recalculating the light curve using each of the three different comparison stars; 
the standard deviation of the results was 0.4 days. Thus, in what follows, we adopt the value 
P TOt = 13.4 ± 0.4 days. Since this result is only based on 2.5 cycles, continued observations are 
warranted to check our estimate. There may be additional errors because of differential rotation 
and the variations in spot positions and intensities on the rotational timescale. The bottom panel of 
Fig. 3 shows the flux as a function of rotational phase during the epoch of quasiperiodic variation. 
A periodogram of the entire data set also shows a peak at 13.4 days, along with a peak at 6.7 days, 
presumably from a time period when star spots occurred on both sides of the star. 

2.3. MAGNUM Observatory 

We observed the transit of UT 2006 Aug 21 with the 2m telescope at the Multi-color Active 
Galactic NUclei Monitoring (MAGNUM) observatory on Haleakala, Hawaii (Kobayashi et al. 1998; 
Yoshii 2002; Yoshii, Kobayashi, & Minezaki 2003). This is the same transit that was observed 
spectroscopically by Winn et al. (2006). We used the multi-color imaging photometer (MIP), 
which allows for simultaneous observation with an optical 1024 2 SITE CCD and an infrared SBRC 
InSb 256 2 detector. In this case, we used only the optical detector, because the infrared detector 
was saturated even in very short exposures. We observed in the V band. Because the MIP field of 
view is only 1'5 x 1^5, and there are no good comparison stars within this small field, we nodded 
repeatedly between the target star and a calibration star, with 1 s exposures. The calibration 
star was HD 190449 (V = 8.12, B — V = 0.79, K0). The average interval between exposures of 
HD 189733 was 46 s. The median FWHM of stellar images was l'/2. 

We reduced the images with the standard MIP pipeline described by Minezaki et al. (2004). 
We then performed aperture photometry on HD 189733 and HD 190449, using an aperture radius 
of 6'.'65 and a sky annulus ranging in radius from 6'.'65 to 9 / /42. To produce a comparison signal, 
the time series for HD 190449 was boxcar-smoothed (with a width of 4 points, or 3 minutes) and 
then linearly interpolated onto the time stamps of the HD 189733 data. Then the HD 189733 time 
series was divided by this comparison signal. A few extreme outlying points were rejected. To 
remove residual systematic errors in the out-of-transit flux determination, we divided by a linear 
function of time that was determined as part of the fitting procedure (see § 4). The final light curve 
is shown in the fifth panel of Fig. 2. 

2.4. Wise Observatory 

We observed the transit of 2006 Sep 5 with the lm telescope at Wise Observatory, in Israel. 
We used a Tektronix 1024 2 back-illuminated CCD detector, giving a pixel scale of / /7 and a field 
of view of llf9 x 11?9. We observed through a Johnson / filter, the reddest optical band available 
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Fig. 2. — Relative photometry of HD 189733 during 6 different transits. The date, telescope, 
and filter are identified on each panel. These data were used only for measurements of the times 
of transit (see § 3). The model, shown as the solid line, is based on the model derived from the 
FLWO z-band data after changing the limb-darkening parameters appropriately. In all cases, the 
residuals (observed— calculated) are plotted beneath the data. 
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Fig. 3. — Relative photometry of HD 189733, from the T10 0.8m APT at Fairborn Observatory. 
The top two panels show data from the 2005 and 2006 observing seasons, respectively. The solid 
symbols show the portion of data that was used in the periodogram analysis. This subset of the 
data is plotted in the bottom panel as a function of the photometric phase. 
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on this camera. The exposure time was 10 s and the telescope was defocused in order to avoid 
saturation. Autoguiding was used to keep the image registration constant throughout the night. 
We also obtained sky flat exposures at sunset and zero-second (bias) exposures at the beginning and 
during the night. We performed the data reduction and photometry using very similar procedures 
to those that were used on the FLWO data (§ 2.1). The bottom panel of Fig. 2 shows the final 
light curve. 

3. Determination of System Parameters 

To estimate the planetary, stellar, and orbital parameters, and the times of transit, we fitted 
a parameterized model to the transit photometry. The model and the fitting method were similar 
to those described in previous TLC papers (see, e.g., Winn et al. 2007), except that in this case we 
accounted for correlated noise, as described below. The model is based on a Keplerian orbit of a star 
(with mass M* and radius R+) and a single planet (M p , R p ) about their center of mass. For most 
of the analysis we assumed that the orbital eccentricity e is zero, because the expected timescale 
of tidal circularization is short in the absence of excitations from other planets (see, e.g., Rasio et 
al. 1996, Trilling et al. 2000, Dobbs-Dixon et al. 2004, Adams & Laughlin 2006). However, in § 4.3, 
we discuss the empirical upper limits on the orbital eccentricity. The orbit has a period P and an 
inclination i relative to the sky plane. We define the coordinate system such that 0° < i < 90°. 
It is often useful to refer to the impact parameter b = a cos i/R* (where a is the semimajor axis) 
rather than the inclination. 

Because one of our goals was to measure the individual transit times, we allowed each transit 
to have an independent value of T c , the transit midpoint, rather than forcing them to be separated 
by exact multiples of the orbital period. Thus, the only effect of P on the model is to determine 
the semimajor axis a for a given value of the total mass. We fixed P = 2.218575 days, the 
value determined by Bouchy et al. (2005) and Hebrard & Lecavelier Des Etangs (2006) from the 
detection of transits in the Hipparcos database. The uncertainty of 0.000003 days was negligible for 
our purposes, although we were able to use the resulting values of T c to produce an independent 
estimate of the period, as described in § 5. 

Neither M* nor M p can be determined from photometry alone. As we have done in previous 
TLC analyses, we fixed at a value that is based on an analysis of the stellar spectrum and other 
observable properties. We then used the scaling relations R p oc M^ 3 and i?* oc M^ 3 to estimate 
the systematic error associated with the uncertainty in M*. In this case we adopted the value 
M* = 0.82 ± 0.03 Mq (Bouchy et al. 2005). The planetary mass M p hardly affects the photometric 
model at all, but for completeness we used the value M p = 1.13 Mj up (Winn et al. 2006). 

To calculate the relative flux as a function of the projected separation of the planet and the 
star, we employed the analytic formulas of Mandel & Agol (2002) to compute the integral of the 
intensity over the unobscured portion of the stellar disk. We assumed the limb darkening law to 
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be quadratic, 

I -f = l-u 1 (l-ti)-u 2 (l-ri 2 , (1) 
n 

where / is the intensity, and [i is the cosine of the angle between the line of sight and the normal to 
the stellar surface. We fixed the limb-darkening coefficients at the values calculated by Claret (2000, 
2004) for observations of a star with the observed spectral properties I§ We also investigated the 
effect of fitting for the limb darkening parameters, as discussed below. In addition, the light curves 
exhibited gradients in the out-of-transit data, probably due to differential extinction between the 
target star and the comparison stars, or some other systematic error. For this reason, each of the 8 
data sets was modeled with two extra parameters: the out-of-transit flux / oot and a time gradient 
a. 



The fitting statistic was 

N f 
3=1 



fj (obs) - fj(calc) lJ 



(2) 



where /j(obs) is the flux observed at time j, crj controls the relative weights of the data points, and 
/j(calc) is the calculated value. In order to derive realistic uncertainties on the parameters, it is 
important for Oj to include not only measurement errors but also any unmodeled systematic effects. 
Of particular importance is the timescale of the systematic effects. Correlated noise effectively 
reduces the number of independent data points and correspondingly increases the uncertainties 
in the model parameters, an issue that Pont et al. (2006) and Gillon et al. (2006) have recently 
brought to attention in the context of transit photometry. 

Our approach to this problem was as follows. First, for each of the 8 transits, we rescaled 
the instrumental uncertainties such that xV-^dof = 1 f° r the best-fitting model. The resulting 
uncertainties are those that are given in Table 2. Second, we followed the procedure of Gillion et 
al. (2006) to decompose the observed noise into "white noise" (that which averages down as 1/y/N, 
where N is the number of data points) and "red noise" (that which does not average down over 
some specified time interval). Specifically, we calculated the standard deviation of the residuals 
(<r) and the standard deviation of the time-averaged residuals (ctjv)- The averaging time was 1 hr 
(a timescale comparable to the transit event), corresponding to a number N of data points that 
depended upon the cadence of observations. Then we solved for the white noise a w and red noise 
a r from the system of equations 

°1 = °"™ + °"r> ( 3 ) 
2 

N = ~W r ' ^' 

Finally, to account approximately for the effective reduction in the number of independent data 
points, we rescaled the aj in Eq. (J2j) by the factor a r /(a w /y/N). 



Specifically, we used the tabulated values for an ATLAS model with T c g = 5000 K, logg = 4.5 (cgs), log [M/H] = 
0.0 and v t = 2.0 km s" 1 . For the z band, Ui = 0.32 and u 2 = 0.27. 
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The results for o~ w , ay, and the rescaling factor for each data set are given in Table 1. These 
results are not very sensitive to the choice of averaging time. Any choice between 15 min and a 
few hours gave similar results. Among the data sets are wide disparities in the degree of red noise, 
ranging over a factor of 10. By far the best data, in the sense of the smallest noise correlations, 
are from the FLWO 1.2m telescope and Keplercam (§ 2.1). For this reason, we decided to estimate 
the system parameters using only the FLWO data, and use the other data sets only to determine 
transit times and as a consistency check on the FLWO results. (Had we included the other data 
sets, their statistical weight would anyways have been much smaller.) 

In short, we used the FLWO data to solve for the two bodies' radii (R± and R p ); the orbital 
inclination (i); and the mid-transit time (T c ), the out-of-transit flux (/ 00 t)) an d a time gradient (a) 
for each of the 2 transits. We then fixed -R*, R p , and i at the best-fitting values, and fitted each of 
the remaining 6 data sets to find T c , f oot , and a. 

To solve for the model parameters and their uncertainties, we used a Markov Chain Monte 
Carlo algorithm (see, e.g., Tegmark 2004). Our jump function was the addition of a Gaussian 
random number to each parameter value. We set the perturbation sizes such that ~20 % of jumps are 
executed. We created 10 independent chains, each with 500,000 points, starting from random initial 
positions, and discarded the first 20% of the points in each chain. The Gelman & Rubin (1992) R 
statistic was within 0.2% of unity for each parameter, a sign of good mixing and convergence. We 
merged the chains and took the median value of each parameter to be our best estimate, and the 
standard deviation as the 1 a uncertainty. 



The results are given in Table 3. Along with the results for the model parameters, we have 
provided results for some useful derived quantities such as the impact parameter b, the radius ratio 
(R p /R+), and the fraction (R p /a) 2 (which gives the fraction of starlight reflected by the planet, 
for an albedo of unity). We also report the calculated values of the full transit duration (the time 
between first and fourth contact, tjy — t\), and the partial transit duration (the time between first 
and second contact, or between third and fourth contact) Q 




The result for the stellar radius is it!* = 0.753 ± 0.025 Rq. The uncertainty is dominated 
by the statistical error of 0.023 Rq (3.1%). The covariance with the uncertainty in the stellar 
mass produces an additional error of 0.009 Rq (1.2%), which we have added in quadrature to the 



4. 



Results 




4.1. Stellar and Planetary Radii 



4 Although the partial transit duration is listed as tn — ti in Table 3, all of the results in Table 1 are based on the 
entire light curves, including both ingress and egress data. Our model assumes tn —ti = tiv — tin. 
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statistical error to arrive at the net uncertainty of 0.025 Rq. We find the planetary radius to 
be R p = 1.156 ± 0.046 i?j U p, where the uncertainty is again dominated by the statistical error of 
0.044 i? Jup . 

To test the robustness of these results, we performed some additional fits. We gauged the 
importance of the choice of limb-darkening law by re-fitting the data under different assumptions. 
When we allowed the limb-darkening coefficients U\ and ui to be free parameters rather than 
holding them fixed, we found i?* = 0.755 Rq and R p = 1.163 -Rj up , well within the 1 a error of 
our original analysis. (The optimized limb-darkening coefficients for the z band were u\ = 0.35 
and U2 = 0.22, as compared to the theoretical values of u\ = 0.32 and U2 = 0.27.) Likewise, the 
results changed by only 0.25 a when we used a linear limb-darkening law (u2 = 0), regardless 
of whether the linear limb-darkening coefficient was fixed or taken to be a free parameter. As 
another check, we fitted each of the 8 transit data sets separately, solving for the parameters 
{R+, R p , i, T c , a, /oot} in each case. Although the statistical power of the FLWO photometry was 
the greatest, as noted previously, all 8 sets of results agreed within their calculated error bars. In 
particular, the unweighted "ensemble averages" of the results from the 6 non-FLWO data sets were 
it!* = 0.760 Rq and R p = 1.148 i?j U p, again in agreement with our original analysis. 

How do these results compare to the previous analyses of HD 189733 by Bouchy et al. (2005) 
and Bakos et al. (2006a)? An important difference is that those authors decided not to determine 
R+ from the transit photometryO Rather, those authors used estimates of R+ based on an analysis 
of other observable properties of the star. Bouchy et al. (2005) used measurements of the star's 
parallax, effective temperature, surface gravity, and metallicity, in comparison with the outputs 
of stellar evolution models, and concluded i?* = 0.76 ± 0.01 Rq. Likewise, Bakos et al. (2006a) 
investigated four different ways of determining the stellar radius, based on broad band colors, 
spectral properties, and model isochrones, and found a stellar radius in the range 0.74-0.79 Rq. 
Ultimately, Bakos et al. (2006a) adopted the value 0.758 ± 0.016 Rq based on a calibration of 
2MASS photometry by Masana et al. (2006). Our determination of R+ based on the transit light 
curve is in agreement with those independent determinations; the mutual agreement constitutes an 
important consistency check on the data and our analysis. 

Our value for the planetary radius agrees with the value R p = 1.154 ± 0.033 Rj up found by 
Bakos et al. (2006a). While the Bakos et al. (2006a) result would appear to be more precise, the 
comparison is somewhat misleading because Bakos et al. (2006a) were not simultaneously fitting 
for R+, as noted above. If we follow their procedure of fixing i?* = 0.758 ±0.016 Rq, then our result 
for the planetary radius becomes more precise: R p = 1.164 ± 0.028 Rj up . Both our result and that 
of Bakos et al. (2006a) disagree with that of Bouchy et al. (2005), who found 1.26 ±0.03 Rj np based 
on a £?-band light curve from the 1.2m telescope at the Observatoire de Haute-Provence (OHP). 



5 While Bakos et al. (2006a) did fit for the stellar radius, finding i?* = 0.68 ±0.02 _Rq, they did not trust the result. 
They quite reasonably suspected that the true error bar was significantly larger than their calculations indicated, 
because of correlated noise in the data. 
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This discrepancy was traced by Bakos et al. (2006a) to systematic errors in the OHP photometry. 
When the B band light curve was recalculated using a greater number of comparison stars, the 
transit depth decreased by 20% and the inferred planetary radius shrank accordingly. 

How does the planetary radius compare to theoretical expectations, given its mass (M p = 1.13± 
0.03 Mj up ; Winn et al. 2006) and its proximity to its parent star (a = 0.035 AU)? Fortney, Marley, 
& Barnes (2006) have recently provided a wide range of theoretical predictions for exoplanet radii. 
Their calculations are for a solar-mass star, but at their suggestion we can rescale the semimajor 
axis to compensate for the lower luminosity of HD 189733. Assuming (L+/L@) = (M+/M & ) s ' 5 , then 
if HD 189733 were orbiting the Sun at a = 0.05 it would receive roughly the same incident flux as it 
does in its actual orbit. The resulting prediction from Fig. 6 of Fortney et al. (2006) is a planetary 
radius between ps1.05-1.12 -Rj up , assuming an age of ~4.5 Gyr and depending on whether or not 
the planet has a massive (25 Mq) core. The core-free prediction for the radius is larger and in 
agreement with the observed value. However, it seems premature to claim that a massive core is 
disfavored, given the uncertainties that enter into the calculations and the uncertainty in the age 
of the system. It does seem safe to say that the radius of HD 189733 does not present a severe 
theoretical problem, unlike the cases of the apparently "bloated" planets HD 209458b, HAT-P-lb, 
and WASP-lb (for recent results on those systems, see Knutson et al. 2006, Bakos et al. 2006b, 
Collier-Cameron et al. 2006, and Charbonneau et al. 2006b). 

The mean density of HD 189733b is p p = 0.91 ± 0.06 g cm" 3 , which is between the densities 
of Saturn (0.6 g cnf 3 ) and Jupiter (1.2 g cm -3 ). The surface gravity of HD 189733b is g = 
21 ± 1 m s~ 2 , which is also intermediate between Saturn (10 m s~ 2 ) and Jupiter (25 m s -2 ). We 
note that whenever M p is measured via the spectroscopic orbit of the star, and R p is measured via 
transit photometry (as is the case here) , then the derived value of g is immune to systematic errors 

2 /3 

in the parameters of the parent star. This is because the fitting degeneracies are M p oc M* and 
R p oc M*^ 3 , and hence g = GM p /R p is independent of M*. Here this fact is of limited interest, 
because the error in R p is dominated by statistical error, but it may be of importance in future 
transit studies I§ 



4.2. Determination of the Transit Ephemerides 

Table 4 gives the 8 transit times measured from our data. We have used these times, along 
with transit times previously measured by Bakos et al. (2006), to calculate a photometric ephemeris 
for this system. Although Bakos et. (2006) reported 15 measured times, we used only 4 of those 
data points in our analyis. We did not include the 10 times that were based on only partial 
observations of the transit. Full transits are greatly preferable, in order to correct (or at least 
assess) systematic errors using the pre-ingress and post-egress data. In addition, we did not include 



6 We thank S. Gaudi for helping us to appreciate this point. 
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the T c measurement based on the OHP 5-band light curve of Bouchy et al. (2005) because of the 
systematic errors noted previously. What remained were 4 data points representing 4 independent 
measurements of the same event. In combination with our 8 data points, we fitted a linear function 
of transit epoch E, 



finding T c (0) = 2453988.80336(23) [HJD] and P = 2.2185733(19) days, where the numbers in 
parentheses indicate the 1 a uncertainty in the final two digits. The fit had xV-^dof = 1-08 
and iVdof = 11- We chose E = to correspond to the most precisely known transit time. Our 
derived period agrees almost exactly with the value 2.2185730(20) days determined by Bakos et 
al. (2006a), and it is also in agreement with the Hipparcos-based values of 2.2185750(30) days 
(Bouchy et al. 2005) and 2.218574±°;™[>) da y s (Hebrard & Lecavelier des Estangs 2006). 

Figure 4 is the O— C (observed minus calculated) diagram for the transit times, according to 
this new ephemeris. The filled symbols represent data points used in the fit. There is not yet any 
pattern in the residuals that would indicate the effect of a perturbing body in the system. The 
unfilled square shows the OHP B-b&nd measurement, which is indeed an outlier. The unfilled circle 
is explained in the next section. 



As mentioned previously, one would expect the orbit of a hot Jupiter such as HD 189733b to 
be very nearly circular, due to tidal effects. Previous results and our results have all shown that 
a circular orbit does indeed provide a satisfactory description of the available data. However, it is 
still interesting to make an empirical determination of the eccentricity, both in the spirit of "what 
one can measure, one should measure" and also because any additional bodies in the system could 
excite the orbital eccentricity. We used a two-step procedure to determine the orbital eccentricity. 

First, we used our revised ephemeris to interpret the secondary-eclipse timing of Deming et 
al. (2006). For a nonzero (but small) orbital eccentricity, the time difference between the midpoint 
of secondary eclipse, T sec , and the time of transit, T tra , may differ from half of the orbital period: 



where u is the argument of pericenter. Deming et al. (2006) measured the midpoint of a secondary 
eclipse at HJD 2453692.62416 ± 0.00067, which is consistent with the e = prediction of our 
ephemeris. The unfilled circle in Fig. 4 represents the secondary eclipse measurement. The timing 
offset of Eq. © is 0.46 ± 1.1 minutes, corresponding to ecosw = 0.00023 ± 0.00054. 

Second, to determine the other component esina; of the eccentricity vector, we used the radial 
velocities presented by Winn et al. (2006). We used only those 60 velocities measured outside of 
transits (i.e., not affected by the Rossiter- McLaughlin effect). We performed an MCMC analysis to 



T C (E) = T c (0) + EP, 



(5) 



4.3. Limits on the Orbital Eccentricity 




(6) 
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Fig. 4. — Transit and secondary-eclipse timing residuals for HD 189733. The calculated times, 
using the ephemeris of Eq. ([5]), have been subtracted from the observed times. The filled symbols 
represent data points used in the fit. The unfilled square is based on 5-band data by Bouchy et 
al. (2005), as re-analyzed by Bakos et al. (2006a). The unfilled circle is the secondary-eclipse time 
measured by Deming et al. (2006). 
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solve for the Keplerian orbital parameters, as well as a possible long-term velocity gradient, using 
the same treatment of the measurement errors that Winn et al. (2006) applied to the entire data set. 
We allowed both e and to to be free parameters, but with an a priori constraint on e cos u to enforce 
compliance with the secondary-eclipse measurement. The result for esina; was —0.007 ±0.011. 

Hence, both components of the eccentricity vector are consistent with zero, and e cos u is about 
20 times more tightly bounded than esinw. Values of e as large as ~0.02 are allowed, but only for 
uj very close to ±90°. 



Thanks to the APT data (§ 2.2), HD 189733 is the first star with a measured rotation period 
that also has a transiting planet. Together with the transit photometry and the observation of 
the Rossiter-McLaughlin effect, this allows for the determination of the angle between the stellar 
rotation axis and the planetary orbit normal, as anticipated by Queloz et al. (2000). This angle is 
worth measuring because any significant misalignment may be an indication of perturbative effects 
during planetary migration, among other reasons (as explained in more detail by Ohta et al. 2005, 
Winn et al. 2005 and Gaudi & Winn 2007). 

The true (three-dimensional) angle if} between the stellar spin axis and the orbital axis is given 
by the formula 



where i is the orbital inclination, is the inclination of the stellar rotation axis, and A is the angle 
between the sky projections of the two axes. (For a diagram of the coordinate system, see Fig. 3 
of Ohta et al. 2005.) The transit photometry determines i with excellent accuracy. Observations 
of the Rossiter-McLaughlin effect have been used to determine A, but they cannot be used to 
determine i* independently. Rather, they are sensitive to vsinz*, the projected rotation rate of 
the star. Given usini* from the Rossiter-McLaughlin observations^, along with i and R± from the 
transit photometry, and the stellar rotation period P ro t, one can determine sini* via the formula 



Hence, all of the angles in Eq. (J7|) are known. Using the values A = — 1?4 ± 1?1, and vsini+ = 
2.97 ± 0.22 km s" 1 from Winn et al. (2006), we found sinz* = 1.04 ± 0.09. By rejecting values 
of sinz* > 1 as unphysical, and propagating the errors through Eq. ([7]), we determined an upper 
bound on the (mis) alignment angle ip of 27° with 95% confidence. Essentially the same result can 
be obtained from the approximation cosip ~ sini*, which is valid because i ~ 90° and A ~ 0° are 
tightly constrained. 



4.4. Three-Dimensional Spin-Orbit Alignment 



cos ip = cos i* cos i + sin i* sin i cos A 



(7) 




(8) 



7 Of course it is also possible to use a more traditional measurement of usinu, from an analysis of the width of 
photospheric absorption lines. 
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This is the first exoplanetary system for which it has been possible to measure ip. The result 
is consistent with zero, but it is not as precise as the result for the projected angle A. How could 
the measurement of ifi be improved? We have already mentioned some caveats relating to the 
measurement of the rotation period^] (§ 2.2), but the current uncertainty in ip is dominated by 
the error in -usini* which is itself dominated by systematic errors arising from the interpretation 
of the transit spectra (see Winn et al. 2006). Specifically, the systematic error arose from the 
"calibration" the Rossiter-McLaughlin effect using simulated spectra, which was needed because the 
transit spectra were observed through an I2 cell and analyzed with an algorithm that is nominally 
designed for measuring Doppler shifts rather than spectral distortions. Further improvement might 
be achieved through a more sophisticated set of simulations or perhaps by re-observing the Rossiter- 
McLaughlin effect without the I2 cell (i.e., using a different technique to account for instrumental 
variations). 



5. Summary 

We have presented photometry of 8 complete transits of the exoplanet HD 189733b, and 
modeled the light curves in order to determine the radii of the star and the planet. Our results 
are consistent with previous results and with theoretical expectations for close-in Jovian planets. 
Stringent limits on the orbital eccentricity follow from the measured transit times, in conjunction 
with a previous detection of the secondary eclipse and with the spectroscopic orbit. We have also 
presented nightly out-of-transit photometry spanning 2 yr that has revealed the stellar rotation 
period. We have used this information, along with a previous analysis of the spectroscopic transit, 
to place an upper bound on the true angle between the stellar rotation axis and the orbital axis. 
With these developments, HD 189733b has become one of the most thoroughly characterized planets 
outside of the Solar system. 

We thank F. Pont and M. Gillon for helpful discussions about correlated noise. We are grateful 
to G. Marcy, P. Butler, S. Vogt, and E. Turner for their help with the Doppler analysis and for 
encouragement. A.R. thanks the MIT UROP office for research funding. 
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Table 1. Characteristics of Transit Data 



Date 


Telescope 


Filter 


Cadence 


White Noise 


Red Noise 


Reweighting 


; factor 


(UT) 






(min) 


o w 


o> 


oy/(<Wv 




2005 Nov 28 


T10 APT 0.8m 


(b + y)/2 


1.44 


0.0045 


0.00250 


4.2 




2006 May 2 


T10 APT 0.8m 


(b + y)/2 


1.44 


0.0024 


0.00080 


2.8 




2006 May 22 


T10 APT 0.8m 


(b + y)/2 


1.44 


0.0024 


0.00090 


3.0 




2006 Jun 11 


T10 APT 0.8m 


(b + y)/2 


1.30 


0.0019 


0.00140 


5.3 




2006 Jul 21 


FLWO 1.2m 


z 


0.23 


0.0023 


0.00035 


3.3 




2006 Aug 21 


MAGNUM 2m 


V 


1.38 


0.0029 


0.00100 


2.9 




2006 Sep 5 


Wise lm 


I 


0.73 


0.0029 


0.00120 


4.6 




2006 Sep 10 


FLWO 1.2m 


z 


0.23 


0.0022 


0.00020 


2.0 





Table 2. Photometry of HD 189733 



HJD Telescope Filter Relative flux Uncertainty 



FLWO 1.2m z 2453937.71893 1.0037 0.0027 

FLWO 1.2m z 2453937.71909 0.9963 0.0027 

FLWO 1.2m z 2453937.71925 1.0007 0.0027 



Note. — The time stamps represent the Heliocentric Julian Date at the 
time of mid-exposure. We intend for this Table to appear in entirety in the 
electronic version of the journal. A portion is shown here to illustrate its 
format. The data are also available from the authors upon request. 
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Table 3. System Parameters of HD 189733 



1— ' O T* O TY1 OT OV 


Value 




(i?,/i? )(M,/O.82 Mq)- 1 ^ 


0.753 


0.023 


( J R P /i?ju P )(M,/0.82 Mq)- 1 ^ 


1.156 


0.044 




0.753 


0.025 


•Rp/ Rjup 


1.156 


0.046 


Rp/ 


0.1575 


0.0017 


(R p /af 


0.000313 


0.000025 


Rja 


0.1124 


0.0034 


i [deg] 


85.76 


0.29 


b 


0.658 


0.027 


tlV - ti [hr] 


1.827 


0.029 


tn — ti [min] 


24.6 


1.9 



Table 4. Mid-transit times of HD 189733b 



Epoch 
E 


Mid-transit time 
[HJD] 


Uncertainty 
[days] 


-129 


2453702.60416 


0.0032 


-59 


2453857.90694 


0.00086 


-50 


2453877.87598 


0.00094 


-41 


2453897.84444 


0.0014 


-23 


2453937.77590 


0.00044 


-9 


2453968.83715 


0.0013 


-2 


2453984.36592 


0.0013 





2453988.80331 


0.00027 



